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Abstract 

We have theoretically studied the full counting statistics of electron transport through a single- 
molecule magnet (SMM) with an arbitrary angle between the applied magnetic field and the SMM's 
easy axis above the sequential tunneling threshold, since the angle 9 cannot be controlled in present- 
day SMM experiments. In the absence of the small transverse anisotropy, when the coupling of the 
SMM with the incident-electrode is stronger than that with the outgoing-electrode, i.e., Tl/Tr ^ 1, 
the maximum peak of shot noise first increases and then decreases with increasing 9 from to 0.57T. 
In particular, the shot noise can reach up to super-Poissonian value from sub-Poissonian value when 
considering the small transverse anisotropy. For Tl/Tr <C 1, the maximum peaks of the shot noise 
and skewness can be reduced from a super-Poissonian to a sub-Poissonian value with increasing 9 
from to 0.57r; the super-Poissonian behavior of the skewness is more sensitive to the small 9 than 
shot noise, which is suppressed when taking into account the small transverse anisotropy. These 
characteristics of shot noise can be qualitatively attributed to the competition between the fast 
and slow transport channels. The predictions regarding of the ^-dependence of high order current 
cumulants are very interesting for a better understanding electron transport through SMM, and 
will allow for experimental tests in the near future. 
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I. INTRODUCTION 



Electronic transport through an individual single-molecule magnet (SMM) has attracted 
intense experimental^ and theoretical^ invention dne to its potential application 
in molecular spintronics devices [141] and classical [15j and quantum information processing 
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171 ]. The prototypal SMM is characterized by a large spin (S > 1/2), easy-axis anisotropy 



which defines the preferred z axis in spin space along which spin is quantized, and transverse 
anisotropies which allow tunneling transitions between the molecular eigenstates of Sf ot . 
Since the transverse anisotropy may lead to mixing of spin eigenstates of Sf ot , the molecular 
eigenstates are not simultaneous eigenstates of S* ot . However, these transverse anisotropy 
terms are very small compared with the easy-axis anisotropy so that they can be taken 
into account by the standard perturbation calculation. On the other hand, the effect of an 
external strong magnetic field on electron transport through the SMM has also been studied, 
in which the easy axis of the SMM is usually assumed to along the direction of the external 
magnetic field B. In the present actual break-j unction and electromigration experiments, 
however, the angle of the external field with respect to the easy axis of the SMM is unknown 
and cannot be controlled[l-4]. If the angle between the easy axis and magnetic field is not 
small, the transverse Zeeman energy may compare with the easy-axis anisotropy energy. 
This implies that the molecular eigenstates are not approximate eigenstates of the spin 
component along any axis, which leads to the failure of the perturbation calculation. In 



very recent single-molecule experiment, Zyazin et a/.[4J found that the angle between B 
and the easy axis of the SMM plays an important role in fitting the theoretical model to 
experimental data. Therefore, it is significant to study the effect of the angle between the 
SMM's easy axis and magnetic field on electron transport in the SMM system. 

Although the present experimental studies focused on the differential conductance or 
average current [l- 4], the full counting statistics (FCS) of electron transport through single- 
molecule magnet or molecular junction has been attracting much theoretical research in- 



terests 



flfl 
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22] owing to its allowing one to identify the internal level structure of 



the transport system 



nig 
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231 ] and to access information of electron correlation that can 



not be contained in the differential conductance and the average current 2411 . For example, 

on 

our previous studies |9|, [lOj have shown that the super- Poissonian noise characteristics of 
electron transport through the SMM can be employed to reveal important information of 
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the internal level structure of the SMM and the left-right asymmetry of the SMM-electrode 
coupling. In addition, the frequency- resolved shot noise spectrum of artificial SMM, e.g., 



a CdTe quantum dot doped with a single S = 5/2 Mn spin, can a 



low one to separately 

extract the hole and Mn spin relaxation times via the Dicke effect [22|. Especially, the FCS 
may provide the full information about the probability distribution P (n, t) of transferring 
n electrons between electrode and SMM during a time interval t. The FCS may be ob- 
tained from the cumulant generating function (CGF) F (x) which related to the probability 
distribution by [25] 

e- F W=^p(n,t)e in3£ , (1) 

n 

where x is the counting field. All cumulants of the current can be obtained from the CGF by 
performing derivatives with respect to the counting field Ck = — (—id x ) k F (x) ■ I n the 
long-time limit, the first three cumulants are directly related to the transport characteristics. 
For example, the first-order cumulant (the peak position of the distribution of transferred- 
electron number) C\ = n gives the average current (J) = eCi/t. The zero- frequency shot 
noise is related to the second-order cumulant (the peak-width of the distribution) S = 
2e 2 C 2 /t = 2e 2 (n 2 — n 2 ^ jt. The third cumulant C 3 = (n — n) 3 characterizes the skewness 
of the distribution. Here, (• • • ) = J2 n (' ' ') ^ ( n ^)- I n g enerai , the shot noise and the 
skewness are represented by the Fano factor F 2 = C 2 /Ci and F 3 = C3/C1, respectively. 

In this work, we consider a more universal model of the SMM and investigate the effect of 
the angle between the easy axis of the SMM and the applied magnetic field and the transverse 
anisotropy on the FCS in SMM. Up to now, the effects of the angle between the easy axis 
and magnetic field, and the transverse anisotropy on the FCS in the present SMM system 
has not been studied to the best of our knowledge. We found that although the threshold 
bias voltage of the sequential tunneling has only a tiny decrease with the increase of the 
angle 6 from to 0.57T, the quantum noise properties of electron transport through SMM is 
not only depend on the left-right asymmetry of the SMM-electrode coupling, but also the 
angle 9 between the easy axis and magnetic field, which can be qualitatively attributed to 
the competition between the fast and the slow transport channels. The paper is organized 
as follows. In Sec. II, we introduce the SMM system and outline the procedure to obtain 
the FCS formalism based on an effective particle-number-resolved quantum master equation 
and the Rayleigh-Schrodinger perturbation theory. The numerical results are discussed in 
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Sec. Ill, where we discuss the effects of an arbitrary angle between the easy axis of the 
SMM and the applied magnetic field, and the second-order transverse anisotropy on the 
super-Poissonian noise, and analyze the occurrence-mechanism of super-Poissonian noise. 
Finally, in Sec. IV we summarize the work. 

II. MODEL AND FORMALISM 

A SMM coupled to two metallic electrodes L (left) and R (right) is described by the 
Hamiltonian H total = H mo i + Hi eads + H T . We assume that the SMM-electrode coupling is 
sufficiently weak so that the electron transport is dominated by sequential tunneling. The 
SMM Hamiltonian is given by 



Here, the first two terms depict the lowest unoccupied molecular orbital (LUMO), n = 
d^d^ + d^di is the number operator of the electron in the molecule, where d) a (d a ) creates 
(annihilates) an electron with spin a and energy Ed (which can be tuned by a gate voltage V g ). 
U is the Coulomb repulsion between two electrons in the LUMO. The third term describes the 
exchange coupling between electron spin in the LUMO and the giant spin, the electronic spin 
operator s = i&o-a') d c i with a = (a x , a y , a z ) being the vector of Pauli matrices. The 

forth and fifth terms are the anisotropy energies of the SMM, where K\ describes the easy- 
axis anisotropy and K 2 the transverse anisotropy. The last term denotes Zeeman splitting, 
where g^B has been absorbed into B. In a general case, since the transverse anisotropy and 
the magnetic field terms do not commute with the easy-axis anisotropy term, S^ ot (= s z + S z ) 
is not conserved and the SMM eigenstates of H mo i are not simultaneous eigenstates of S^ ot . 

— * 

On the other hand, if the external magnetic field B is applied along the easy-axis, in the 
absence of the transverse anisotropy the eigenvalue m of S% ot is a good quantum number, 
hence allowing us to numerically diagonalize the molecular Hamiltonian H mo i in the basis 
represented by the eigenvalue m of Sf ot and the corresponding occupation number n of the 
LUMO level, i.e., (|0,m) , \i) \m - \) , |t) \m+ |> , |2,m)}, where m G [-S,S]. 

The relaxation in the electrodes is assumed to be sufficiently fast so that their electron 
distributions can be described by equilibrium Fermi functions. The electrodes are modeled 
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Hmoi = (e d - eV g )h + —n(n — 1) — Js-S 





as non-interacting Fermi gases and the corresponding Hamiltonian 

H Leads = ^ ^akcr ( ^ Q ,i {0 -Q'akg-; (3) 
akcr 

where a^ k(T (a a k CT ) creates (annihilates) an electron with energy £ a kcr, momentum k, and 
spin a in a (a = L, R) electrode. The tunneling between the LUMO and the electrodes is 
described by 

H T = J2(tA k J° + H.c). (4) 

akcr 

In sequential tunneling regime, the transitions are well described by quantum master 
equation of a reduced density matrix spanned by the eigenstates of the SMM. The detailed 
derivation of the FCS based on the particle-number-resolved quantum master equation can 
be found in Refs. 26H28| . and here, we only give the main results. Under the second order 
Born approximation and Markovian approximation, the particle-number-resolved quantum 
master equation for the reduced density matrix is given by 

fF> it) = -iCp^ it) - ifcpW (t) , (5) 

with 

K P (n) (t) = J2 [44"V B) (*) + p in) (*) 4 +) 4 

-4> (n_1) (*) 4 - 4p (n+1) (*) + ff-c (6) 

where = J2 a =L,R A< $i A <w = ^a n t d ^ n t = fa,n~ = 1 - f a (f a is the Fermi 
function of the electrode a), and T a=LR = 2irg a=L>R |t Q =L,_R| 2 - Liouvillian superoperator C 
is defined as C (• • • ) = [i/moj, (•••)]> an d 9a=L,R are the density of states of the metallic 
electrodes, p^ (t) describes the reduced density matrix of the SMM conditioned by the 
electron numbers tunneling through the right junction up to time t. Throughout this work, 
we set e = ft = 1. Here, the validity of the Markovian approximation deserves some 
discussions. For the case of sequential tunneling, the Markovian approximation is valid when 
the system conductance is small compared to the quantum conduct ance[29j, i.e., I/V <C 
e 2 / (2nh) = 1/ (27r), here we have utilized e = K — 1. In the present SMM system, the 
value of I/V is of the order of 10~ 3 < 1/ {2ir). This means that the typical time between 
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two tunneling events is r = e/I ^> 2nh/ (eV) = Ti ea d s , i.e., the SMM dynamics is indeed 
much slower than the decay of lead correlations, thus, the Markovian approximation is 



well justified[29|. The CGF connects with the particle-number-resolved density matrix by 
defining S (x, t) = Y^n P (t)e mx . Evidently, we have e~ F<yX ^ —Ti[S (x, t)), where the trace 
is over the eigenstates of the SMM. Since Eq. (j5J) has the following form 



P 



(n) = Ap (n) + Cp (n+1) + D p (n-l) ^ ^ 



then S (x, t) satisfies 

S = AS + e~ ix CS + e ix DS = L X S, (8) 

where the specific form of L x can be obtained by performing a discrete Fourier transformation 
to the matrix element of Eq. ()5]). Here, the master equation contains off-diagonal matrix 
elements p^ M , which corresponds to superpositions between molecular eigenstates \m) and 
\n). In fact, since the presence of noncommuting Zeeman and transverse anisotropy terms in 
the SMM Hamiltonian, any two eigenstates differ in the spin expectation value (Stoi), which 
leads to different long-range (dipole) fields. Thus the unavoidable interactions between the 
SMM and many degrees of freedom in the environment (e.g., electron spins) lead to rapid 



decay of superpositions of these eigenstates and thus of pfJ^ M |8|, [29|, [30j. As a result, in the 
following calculation the off-diagonal matrix elements can be neglected, and it is sufficient 
to consider the diagonal components of p SMM . 

In the low frequency limit, the counting time (i.e., the time of measurement) is much 



longer than the time of tunneling through the SMM. In this case, F (x) is given by 25l. l3lM33| 



F(x) = -^i(x)t, (9) 

where Ai (x) is the eigenvalue of L x which goes to zero for % — > 0. According to the definition 
of the cumulants one can express Ai (x) as 

^w = £t¥- < 10 > 

fe=i 

Low order cumulants can be calculated by the Rayleigh-Schrodinger perturbation theory in 
the counting parameter x- in order to calculate the first three current cumulants we expand 
L x to third order in x 

L x = L + L 1 x+^L 2X 2 + ^L 3 x 3 + --- ■ (11) 



Along the lines of Ref. 3l|, we define the two projectors V = V 2 = |0)) ((0| and Q = 
Q 2 = 1 - V, obeying the relations VL = L V = and QL = L Q =L . Here, |0)) 
being the steady state p stat is the right eigenvectors of L , i.e., L |0)) = 0, and ((5| = 1 
is the corresponding left eigenvectors. In view of Lq is regular, we can also introduce the 
pseudoinverse according to R — QLq X Q, which is well-defined because the inversion is 
performed only in the subspace spanned by Q. After a careful calculation, Ai (x) is given 
by 



Ai(x) = ((0| J Li|0))x 

+ ^K<0|^|0))-2«6|L 1J RL 1 |0))] x 2 

+ ^ [«0 1 L 3 1 0)) - 3 «0 1 (L 2 RL 1 + L X RL 2 ) \ 0)) 

-6 ((0\ UR^RUP - UR) U 1 0>>] X 3 + --- ■ (12) 



From Eqs. (TTOT) and ffl2l) we can identify the first three current cumulants: 

C 1 A=«0|L 1 |0»/i J (13) 
C 2 /t = [((0 I L 2 1 0)) - 2 «0 | LiRLi \ 0))] /* 2 , (14) 

C 3 /t = [((0 I L 3 1 0)) - 3 «0 | [L 2 RLx + L X RL 2 ) \ 0)) , 
-6 <<0 | L X R {RL X P - L X R) L x \ 0))] /i 3 . (15) 

The three equations above are the starting point of the calculation in following. 

III. NUMERICAL RESULTS AND DISCUSSION 

We now study the effects of the angle of the external field with respect to the easy axis 
of the SMM and the transverse anisotropy on the FCS of electronic transport through the 
SMM weakly coupled to two metallic electrodes. We assume the bias voltage (Vf, = Hl~ /Ar) 
is symmetrically entirely dropped at the SMM-electrode tunnel junctions, which implies that 



the levels of the SMM are independent o: 
not symmetric. Since our previous work 



the applied bias voltage even if the couplings are 



lOf ] has studied the effect of Coulomb interaction 



U on FCS in the SMM in the absence of the transverse magnetic fields and transverse 



anisotropy, we here take a fixed value of U. The parameters of the SMM are chosen as 



B 



80r, where V is the typical 



S = 2, e d = 200r, U = 100r, J = lOOr, K x = 40r and 
tunneling rate of electrons between the SMM and the electrode. In the present work, we 
only study the transport above the sequential tunneling threshold, i.e., > 2e se , where 
e se is the energy difference between the ground state with charge N and the first excited 
states N — 1[34J- In this regime, the inelastic sequential tunneling process is dominant, thus 
electrons have sufficient energy to overcome the Coulomb blockade and tunnel sequentially 
through the SMM. Here, it should be noted that since in the Coulomb blockade regime the 
current is exponentially suppressed and the electron transport is dominated by cotunneling, 
when taking into account cotunneling the normalized second and third moments will deviate 



351 ] . In this paper, we put emphasis 



from the results obtained by only sequential tunneling 
on the effects of the angle between B and the easy-axis of the SMM, and the transverse 
anisotropy on super-Poissonian noise for large left-right asymmetry of the SMM-electrode 
coupling. 

Since the transverse anisotropy and the transverse component of the applied magnetic 
field can lead to mixing of spin eigenstates of S^ ot , the transitions, which are inhibited due 
to spin selection rules in the absence of a transverse field and transverse anisotropy, may 
occur. In order to show explicitly the effect of the angle between B and the easy axis on 
electron transport, we first neglect the small transverse anisotropy. In this case, the applied 
magnetic field B may be assumed to lie in the xz plane because of the rotational symmetry 
of H mo i. Moreover, it is helpful to analyze the selection rules for the occurrence of the 
sequential tunneling. In the absence of the transverse fields and the transverse anisotropy, 
the eigenvalue m of S^ ot is a good quantum number and the sequential tunneling requires 
a change of the electron number by An = ±1, and the magnetic quantum number by 
Am = ±1/2. But for the present case, the only selection rule An = ±1 is still valid, which 
means that arbitrary two states satisfying An = ±1 can do sequential tunneling. For local 
large spin S, there are 25* + 1 empty molecular states with n = in the LUMO, 2 (25* + 1) 
singly-occupied molecular states with n = 1 and 25* + 1 doubly-occupied states with n = 2. 
Therefore, there are 4(2S'+1) 2 transitions, namely, 2 (2S + l) 2 between molecular states 
with n = and n = 1, and 2 (2S + l) 2 between molecular states with n = 1 and n = 2, 
which leads to a much more complex electron transport channels than the case without 
the transverse anisotropy and transverse field. For this reason, we focus on studying the 



S 



dependence of the maximum noise values on the angle 9 in sequential tunneling regime. 

When the coupling of the SMM with the left electrode is stronger than that with the 
right electrode, i.e., Y L /Y R ^> 1, here we choose Y L /Y R = 10. Figures l(a)-(c) show the 
average current, shot noise and skewness as a function of the bias voltage for 9 = 0, O.Itt, 
0.27T, 0.37T, 0.47T, 0.57T. Since the FCS for 9 = 9q has the same bias-voltage-dependence 
as that for 9 = tt — 9q, which arises from the symmetry of the SMM Hamiltonian, we 
restrict our discussion to the case of 9 <G [0,0. 57r]. With increasing 9, the corresponding 
sequential tunneling threshold bias voltage has a tiny decrease and reach their minimums 
when 9 increases to 0.57T, see Fig. 1(a); but quantum noise obviously depends on the angle 
9, see Fig. 1(b) and (c). The maximum peak of shot noise firstly increases and then 
decreases with increasing the angle from to 0.57T. This characteristics of the shot noise 



can be understood with t 
transport channels 
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re help of the dynamic competition between effective fast and slow 



,128| 
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36M38I]. The molecular channel current is given by 



as 



I\n,i) — >\n-l,j) 



C\n-l,j),\n,i)^Rn R \€\n,i) ~ e |n-lJ) — Hr) P\n,i), (16) 
-l,i> — >\n,i) 

— C\ n -i t i) t \ nt j)YRn R \€\ n ,i) — e \n-l,j) — Hr) P\n-l,j)- (17) 



Here C\ n -.xj),\n,i) = \(n — 1, j\ d a \n, i)\ 2 is a constant which related to the two molecular 
states but independent of the bias voltage, where \n,i) {i = 01, 02, 03, 04, 05, for S = 2) 
denote the eigenstates of the molecule with n electrons tunneling into the molecule, which are 
arranged in an ascending order of their eigenvalues e\ n ^). P\ n ,i) is the occupied probability of 
the state \n,i). Since the maximum value of shot noise appears at a large bias voltage, the 
Fermi function f R [e\ n ^ — t\ n -i,j) — Hr) changes very slowly with increasing bias voltage, 

i- e - ; fR { e \n,i) — e \n-i,j) — ^r) — 0- Thus the molecular channel currents I\ n ^ y\ n -i,j) are 

mainly determined by the probability distribution P\ n ,i), and I\ n -i,j) >\n,i) — 0. In the 

presence of the transverse field and the transverse anisotropy, since the transitions between 
the molecular eigenstates are not restricted by the selection rule Am = ±1/2, the possible 
transport channels are 4 (2S + l) 2 , for example, there are a hundred transport channels for 
5 = 2. Therefore, it is unpractical to give all the channel currents. In order to give a 
qualitative explanation for the effect of the angle 9 on the shot noise, we plot the occupied 
probability of the five main molecular eigenstates as a function of bias voltage Vj, for 9 = 0, 



0.17T, 0.27T, 0.37T, 0.47T, 0.57T in Fig. 2. For the case of K 2 = 0, the increase (or decrease) 
of the probability of the molecular eigenstate with high occupancy is always accompanied 
by the decrease (or increase) of the probability of the molecular eigenstates with the low 
occupancy. It is important that the active competition between the fast and slow transport 
channels depends on the angle 9. The competition between the probability of the five main 
molecular eigenstates for 9 = 0.2n, 0.3n and OAn, as shown in Fig. 2, is stronger than that 
for 9 = and 0.57T. Thus, the corresponding transport channel currents can form the so- 
called effective fast-and-slow transport channels, which leads to the maximum value of shot 
noise for 9 = 0.2n, 0.3n and OAn are larger than that for 9 = and 0.57T. In addition, it is 
interesting that some certain angles 9 (e.g., 9 = OAn, OAn, 0.5n) may decrease the maximum 
super-Poissonian value of the skewness F 3 > 1 to sub-Poissonian value of F 3 < 1 although 
the angle (e.g., 9 = 0.2n) can also increase the maximum skewness value, see Fig. 1(c). 
In the simultaneity presence of the transverse field and the small transverse anisotropy, the 
active competition is further strengthened, so that the maximum shot noise value, which is 
sub-Poissonian value for 9 = 0, may be enhanced to super-Poissonian value for some certain 
angles, e.g., 9 = 0.2n, 0.3n and 0.47T, see Fig. 1(e). 

Compared to the case of T L /T R ^> 1, for T L /T R <C 1 the maximum shot noise and the 
skewness peaks are suppressed with increasing the angle 9. Fig. 3(a)-(c) show the average 
current, shot noise and skewness as a function of the bias voltage for 9 = 0, OAn, 0.2n, 0.3n, 
OAn, 0.5n at T L /T R = 0.1. In this situation, the sequential tunneling threshold e se has the 
same characteristics as the case of Tl/T r ^> 1. Apart from this feature, another important 
finding is that with increasing the angle 9 from to 0.5-7T, the maximum peaks of the shot 
noise and skewness can be reduced from super-Poissonian to sub-Poissonian value in the 
absence of the transverse anisotropy, see Fig. 3(b) and (c). Especially for the skewness, its 
super-Poissonian behavior seems more sensitive to the small 9 than shot noise, see Fig. 3(c). 
The shot noise characteristics can also be understood in terms of the so-called fast and slow 
transport channels mechanism. Fig. 4 shows the occupancy probability of the five main 
singly-occupied molecular eigenstates as a function of bias voltage Vb for various values of 
the angle, which determine corresponding transport channel currents. With increasing the 
angle 9 from to 0.57T, the competition between the probabilities of the eigenstates with 
high occupancy and the eigenstates with the low occupancy is gradually weakened, see Fig. 
4. This means that the active competition between the fast-and-slow channel currents is 
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gradually suppressed, thus leading to the maximum super-Poissonian value of shot noise 
is reduced even to sub-Poissonian value. Moreover, the small transverse anisotropy for 
^l/Fr <C 1 also suppresses the maximum peaks of the shot noise and the skewness, and 
thus effaces the sensitivity of the maximum peaks of the shot noise and the skewness to the 
small 9, see Fig. 3(e) and (f). 



IV. CONCLUSIONS 



We have studied the FCS of electron transport through a SMM with an arbitrary angle 
between the external magnetic field and the easy axis of the SMM above the sequential 
tunneling threshold. Since the presence of the transverse field and the transverse anisotropy 
leads to mixing of the spin eigenstates of S^ ot , the spin selection rule Am = ±1/2 for 
sequential tunneling processes are no longer applied to our model, as a result, there are 
4 (2S + l) 2 transport channels participating in the electron transport. Therefore, the angle 
has a complex impact on the FCS. To facilitate the discussion of the origin of the shot 
noise, we put special emphasis on the dependence of the maximum noise on the angle of 
external magnetic field for strong asymmetric coupling to the two electrodes. For the case 
of Fl/Tr ^> 1, the maximum peak of the shot noise firstly increase and then decrease with 
increasing 9 from to 0.57T. In particular, the shot noise is further enhanced and even 
reaches super-Poissonian value when considering the small transverse anisotropy. For the 
case of Tl/Tr <C 1, the maximum peaks of the shot noise and skewness can be reduced 
from super-Poissonian to sub-Poissonian value with increasing the angle 9 from to 0.57T. 
Especially for the skewness, its super-Poissonian behavior seems more sensitive to the small 
angle 9 than shot noise, but this feature is suppressed when taking into account the small 
transverse anisotropy. These characteristics of shot noise can be understood as a result 
of the active competition between the fast and slow transport channels. The predictions 
regarding the high order current cumulants are very interesting for better understanding 
electron transport through individual single-molecule magnet, and the ^-dependence of the 
FCS can be helpful in understanding the experimental results because the angle 9 is difficult 
to be controlled experimentally. 
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FIG. 1: (Colour online) The average currant, shot noise and skewness versus bias voltage for 
different angles of external magnetic field with Fl/Tr = 10. (a), (b) and (c) for Ki = 0, (d), (e) 
and (f) for K 2 = O.LKi. The molecular parameters: S = 2, e d = 200r, U = WOT, J = 100T, 



K\ = 4or, 



B 



80r and k B T = WT. 
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FIG. 2: (Colour online) The probability distribution of molecular eigenstates versus bias voltage for 
different angles of external magnetic field with Tl/Tr = 10 and K2 = 0. The molecular parameters 
are the same as in Fig. 1. 



15 




FIG. 3: (Colour online) The average currant, shot noise and skewness versus bias voltage 
different angles of external magnetic field with Fl/Fr = 0.1. (a), (b) and (c) for K2 = 0, (d), 
and (f) for K2 = O.LKi. The molecular parameters are the same as in Fig. 1. 
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FIG. 4: (Colour online) The probability distribution of molecular eigenstates versus bias voltage 
for different angles of external magnetic field with Tl/Tr = 0.1 and K2 = 0. The molecular 
parameters are the same as in Fig. 1. 
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